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ABSTRACT 



An existing computer simulation of a hydraulic submarine steering 
control system was analyzed to determine the source of its unrealistic 
behavior, which included verifying mathematical equations and numerical 
parameters, duplicating the simulation on another computer, and investi- 
gating the validity of the assumptions made in its development. Rudder 
turning rate and load pressure drop oscillations which were thought to 
be unrealistic were caused by the neglect of significant rudder moments 
related to the control surface turning rate. Knowledge of actual system 
performance was found to be incomplete and further research was suggested. 
Potentially beneficial modifications in the servovalve and actuator 
simulators were identified. Continuing work in the field of unsteady 
rudder behavior was suggested in order to improve simulator and actual 
system designs. 
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I. THE NATURE OF THE PROBLEM 



A. BACKGROUND INFORMATION 

In 1971, the Naval Ship Research and Development Center (NSRDC), in 
Annapolis, Maryland, identified a potential stability problem in an 
electrohydraulic servosystem which was being developed. The control 
system was being designed to control the movement of a submarine's rudder 
in accordance with commands from the helmsman. The NSRDC developed an 
analog computer simulation (hereafter referred to as the "baseline simu- 
lation") which incorporated many of the characteristics of an actual 
hydraulic control system, including friction, hydraulic and mechanical 
characteristics, servovalve behavior, and interdependence of elements 
necessary to create a simulation of the actual system. The simulation 
was based on equations developed to model an existing system used in 
the USS PARGO (SSN 650) . Provisions were made so that the simulation 
could be adjusted by changing values of the various components of the 
model, in order to conform to possible changes of the parameters of 
the system under development, such as changes in dimensions, materials, 
and operating conditions. The provision for adjustment of the simulator 
parameters also allowed simulation of other existing hydraulic systems 
of this type. All parameters which were assumed to be of significant 
impact on the performance of the system were incorporated into the 
simulation. 

When the simulator parameters were arranged to model the performance 
of the system under development, the model did not conform to expected 
system behavior; oscillations in pressures and control surface turning 
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rates were observed which had not been observed in similar systems in 
use. Further investigations failed to reveal any obvious errors in the 
simulation, or deficiencies in the parameters which were believed to be 
of significant impact on the performance of the system. 

This thesis is a result of a project initiated in support of the 
simulation attempt at the NSRDC. The purpose of this project was to 
determine the reasons for the observed behavior of the hydraulic system 
model, and to develop appropriate modifications and refinements of the 
baseline simulation, in order to more closely model actual system 
performance. 

B. SYSTEM UNDER STUDY 

1 . Introduction 

The following paragraphs describe the general characteristics of 
the idealized system that the NSRDC used in the development of its simula- 
tion. The description does not completely describe any real system. 

Rather, it is a model which was designed to include all of the steering 
system characteristics that were expected to be of significant impact on 
the performance of the system. 

2. General System Descriptions 

The steering system (Figure 1) included the servovalve, connecting 
tubing, actuating ram, mechanical linkages, rudder, and minor components. 
The operation of the system is described below. (Numbers in parentheses 
refer to corresponding labeled parts on Figure 1.) 

A desired change of rudder angle is initiated by moving the control 
wheel (1) by hand, which generates an electrical signal having a voltage 
proportional to the degree of movement. This voltage representing the 
desired control surface angle is then compared with a voltage proportional 
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to the actual position of the hydraulic ram (6) . An error signal is 
generated which approximates the difference between the desired and 
actual rudder positions (2) . This signal actuates an electric torque 
motor (3) in the servovalve, which moves a small flapper (4) located in 
the narrow space between two nozzles, causing an increase in the pilot 
pressure on one nozzle, and a decrease on the other. This difference 
in pressure causes the spool valve (5) to move, thereby permitting 
pressure and subsequent fluid flow to be applied to the hydraulic ram 
(6) , resulting in a correcting movement of the control surface (7) through 
the mechanical linkage. The value of the error signal is modified by 
movement of the hydraulic ram, and the entire process is continuously 
repeated until the error signal is reduced to zero. As the error signal 
decreases, the servovalve closes, control surface movement ceases, and 
equilibrium is established at the new position. 

4 

a. Component Descriptions 

(1) Servovalve . The servovalve (Figure 2 ) used in a submarine 
steering control system is a two-stage, force feedback model, incorporating 
a torque motor actuated by an electrical error signal; a flapper-type 
pilot valve, the flapper itself also serving as the armature of the torque 
motor; and a spool valve, actuated by the fluid flow caused by the movement 
of the pilot flapper valve. Movement of the spool results in the direction 
of fluid flow to the actuating ram. 

(2) Actuator . The actuator is a linear, two-way hydraulic 
ram. The NSRDC assumed in modeling the actuator that it was of compensated 
construction; that is, it was designed with equal effective areas on both 
sides of the piston. O-rings are incorporated into the design of the actua- 
tor, eliminating essentially all leakage from one side of the piston to the 
other. 
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CCD 
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FIGURE 2 
SERVOVALVE 



(3) Mechanical Linkage . The linkage connecting the hydraulic 



ram to the rudder consists of a long shaft, support bearings, a packing 
at the penetration of the pressure hull, a crosshead at the end of the 
shaft opposite the ram, and a rod connecting the crosshead to the rudder 
tiller. The crosshead provides a supported, pinned joint from which the 
connecting rod originates. It was assumed to be constrained from lateral 
movement by an open-ended sleeve which was flooded with sea water. 

(4) Rudder . The rudder is composed of two portions, an upper 
rudder and a lower rudder, interconnected to the tiller by a single shaft, 
and located on the top and bottom of the ship's hull, respectively. The 
rudder is symmetric along the chord length, and is designed for operation 
within the range of left 35 degrees rudder to right 35 degrees. 

(5) Other Components . Other components include piping; an 

electric transmitter associated with the control wheel; an error feedback 

* 

servo associated with the actuator; a hydraulic pressure supply which is 
assumed for design purposes to be constant in pressure, and not subject 
to influence by steering system conditions; and a petroleum-base hydraulic 
fluid. 

C. NSRDC BASELINE SIMULATION 

The following sections describe the work done at the NSRDC in order 
to simulate the submarine steering control system. The governing equations 
for each of the components are stated as they were derived by the NSRDC, 
along with any assumptions that were explicitly stated. A schematic circuit 
diagram of the analog computer arrangement is presented, and comments from 
the NSRDC about the nature of the problem and suspected origins of the 
oscillation are stated. 
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1 . Governing Equations and Assumptions 



The governing equations that were provided for each of the compo- 
nents are listed verbatim in this section. Only the governing equations, 
assumptions, and values for constants were communicated to the writer. 
Consequently, the procedure by which the expressions were derived, whether 
through experimentation, experience- justified simplification of compli- 
cated transfer functions, or some other way, is unknown. See Appendix A 
for explanations of the notation used, 
a. Servo Valve 

The servo valve equations include expressions for the valve 
flapper displacement, spool valve displacement, nonlinearity compensations, 
and fluid flow rates through the valve. 

(1) Valve Flapper Displacement . The valve flapper displace- 
ment, x, is expressed as 



x 



sat i 0C 1 (9 o -e)-K 3 t) 
2 



+ .003 
-.003 



( 1 ) 



where is the amplifier gain coefficient, K ^ is the flapper centering 
spring rate, 6 q is the desired position of the rudder, 9 is the actual 
position of the rudder, is the wand spring rate, and z is the present 
valve spool position. No explicit assumptions were given with regard to 
this equation. 

(2) Valve Spool Displacement . The valve spool displacement 
is obtained through the variation with time of the valve spool velocity z: 



z = 



( 2 ) 



0.5 z>. 5 and SGN(z) = SGN(x) 

K^x otherwise 

where is the pilot-stage valve flow coefficient. This expression is 

/ 
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best explained by stating that the spool velocity is proportional to the 
flapper valve displacement, until the spool is fully displaced against 
the valve body, where the movement must of course cease. The velocity 
again becomes non-zero only if the direction of movement is reversed, and 
the valve is allowed to travel away from the hard restraining block. 

Nonlinear effects in the valve, such as varying flow gain 
(i.e., flow is not proportional to spool displacement, due to the shape 
of the ports used in the spool valve), spool friction, deadband, and torque 
motor hysteresis are included in the next equations. 

Dead band is included in the system through the implementa- 
tion of the following equations: 

( z - db where z >db, z>0 

z + db where |zl >db, z<0 (3) 

0 if -db<z<db 

where z' is the valve spool displacement after adjustment for spool dead 
zone. These expressions give the effective spool displacement, i.e., 
adjust the spool displacement to the value it would have for the given 
flow if no dead zone were constructed in the valve. 

Multiple flow gain is accommodated by combining z' with 

z'l z'l : 

z" = .5z' + .Sz'l z'l . (4) 

where z" is the valve spool displacement after spool dead zone, and after 
multiple flow gain. 

Compensation for other nonlinearities such as spool 
friction and torque motor hysteresis may have been accomplished through 
adjustment of the coefficients in the equations for spool dead band, and 
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multiple flow gain. No additional expressions for these nonlinearities 
were given. 

(3) Flow Equation . The flow rate in the spool is related to 
the corrected spool displacement z", the load pressure drop P^, and the 
supply pressure P g by the following expression: 



Q = z"K 5 sgn(P s -sgn(z)P L ) y P s - sgn(z)P L (5) 

where is the main-stage valve flow coefficient. No explicit assumptions 
were listed for this equation, 
b. Actuating Ram 

The position of the actuating ram is described in two ways; by 
the position of the rudder itself, and by the integral of the flow rate 
which has passed through the valve. 

(1) Ram Position in Terms of Rudder Position . The ram position 
y is related to the rudder position, with correction for mechanical resilience 
in the linkage, by 



y = LsinG 




( 6 ) 



where L is the length of the tiller, 0 is the rudder angle in degrees, 

is the mechanical resilience coefficient of the system, and A is the actuator 
piston area. 



(2) Ram Velocity Expressed in Terms of Pressure Change and Flow . 
The NSRDC determined that the ram velocity, flow rate, and load pressure 
have the following relationship to one another: 

Q = -Ay + k 6 p l (7) 



where is the system hydraulic resiliance, and P^ 
of load pressure drop. 



is the rate of change 
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In the development of these equations it was assumed that 
there was no leakage from the control valve to the actuating cylinder at 
line connections, seals, and so forth, and that there was no internal 
leakage in the actuating cylinder itself. Also, it was assumed that a 
compensated actuator was used; that is, an actuator was used in which both 
areas, on opposite sides of the piston, were equal, 
c. Forces and Torques on the System 

All forces on the system were expressed in terms of the equiva- 
lent torques which were placed on the rudder structure. A moment balance 
can then be written: 

M = Tp + T h + T F C8) 

where I is the polar moment of inertia of the rudder mass and the included 
water mass, 0 is the angular acceleration of the rudder, Tp is the torque 
developed as a result of the pressure drop across. the actuator, Tj^ is the 
hydrodynamic torque, and Tp is the friction torque. Each torque will now 
be discussed in more detail. 

(1) Friction Torque . Tp is the most complicated of the torques 
to describe. It is the combined friction torques from the shaft seals, the 
resultant of the normal force on the rudder bearing surface, and any other 
sources of sliding (coulomb) friction present in the system. Stationary 
objects in contact with each other stick together until subjected to forces 
larger than the friction force present between the objects. Additionally, 
the level of applied force necessary to sustain movement at very low speeds 
usually exceeds the force needed at higher velocities. In order to model 
.these system characteristics, the NSRDC developed the following expressions 
for Tp: 
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T„ = 



<VV if |t h *t p |<t s 

) -T s sgn(T H .T p ) if |VTjj >Tj) 
T c sgn(6) le|>e 



0 < e 



( 9 ) 






where e is theoretically zero, but in practice is set to approximately 
0.2 percent of the maximum angular velocity of the rudder. = kTg, 
where k is slightly less than one, typically 0.95, and Tg is expressed by 



T S = c d r|*w + y(K 9 0) 2 + K n e+A-P L ) 2 



+ K 



10 



( 10 ) 



where Cp is the coefficient of friction of the stock bearing, W is the 
weight of the rudder, R is the rudder stock radius, K 9 is the lift co- 
efficient of the rudder, K^ is the drag coefficient of the rudder, and 
Kjq is the seal friction coefficient. It should be noted that the NSRDC 
equation for friction torque did not include any terms for damping pro- 
portional to angular velocity of the rudder surface. 

(2) Hydrodynamic Torque . is the torque generated by the 
rudder when it is placed at some angle of attack other than zero, due to 
the fact that the center of pressure on a rudder does not coincide with 
the turning axis of the structure. The hydrodynamic torque is expressed 
as 

2 

T H = (y ( K 7 9 ' lt 8 e|91 5 tn) 



where V is the ship velocity, is the maximum ship velocity, and Ky and 

K are hydrodynamic torque load coefficients. By adjusting the values of 
8 

Ky and Kg, the hydrodynamic torque variation with rudder position can be 
adjusted to approximate almost any rudder design. If Ky and Kg are 
redefined, the above expression can be written in nondimensional form: 
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HYDRODYNAMIC TORQUE Jh 




FIGURE 3 

HYDRODYNAMIC TORQUE BEHAVIOR 
FOR VARIOUS VALUES OF K y AND Kg 



22 



( 12 ) 



H 



Hmax 





0 _ 

0 

max 



where is the maximum hydrodynamic torque corresponding to Vq, and 

®max t ^ le max i mum rudder angle. See Figure 3 for curves of 
for varying values of K^ and Kg. It can be seen that a variety of rudder 
torque characteristics can be simulated by adjustment of K 7 and Kg. 

(3) Pressure Torque . The pressure torque Tp is the torque 
that is placed on the stock of the rudder by the hydraulic ram, through 
the connecting linkage. It can be expressed as 



Tp = LAP l cos 0 , (13) 

where L is the tiller length, A is the area of the actuator piston, and 
P^ is the load pressure drop. 

2. Analog Computer Simulation 

4 

a. Hardware and Parameters 

Figure 4 is a schematic diagram of the analog computer program 
which was designed to model the governing equations. Table I gives a list 
of typical system dimensions and other parameters which were used to set 
simulation values. Table II gives a list of the potentiometer settings 
for the values listed in Table I, and also lists the expressions which 
were used to develop the settings. 

3. Discussion of the Baseline Simulation 
a. Performance Curves 

Examination of sample recordings obtained from the NSRDC 
showed that oscillations in velocity and pressure were present in the 
simulation outputs. The NSRDC indicated that the oscillation appeared 
to be a function of the inertia of the steering linkage, ram, and rudder, 
and of the system mechanical and hydraulic resiliences. Variation of 
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SERVO VALVE 
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FIGURE 4 

BASELINE ANALOG COMPUTER SCHEMATIC DIAGRAM 




Table I 



Submarine Steering System Parameters Used in the Baseline Simulation 



Parameter 


Value 


Parameter 


Value 


A 


80 in 2 


L 


30 in 


C d 


.05 


P 

r max 


6,000 psi 


db 


.05 


p 

Lmax 


3,000 psi 


I 


9100 in-lbf-sec 2 


P S 


3,000 psi 




.6 


Qmax 


278.86 in^/sec 


K 2 


110 


R 


18 in 


K 3 


5.0 


T Fmax 


1 . 1661X10 6 in-lbf 


K 4 


375 


T 

Hmax 


3.93X10 6 in-lbf 


% 


7.20 


T 

Pm ax 


7.2X10 6 in-lbf 


K 6 


0.015 


W 


25,000 lbm 


K 7 


1.5 


x 

max 


.002 in. 


K 8 


2.5 


y 

'max 


17 in. 


K 9 


30,000 


z 

max 


.5 in. 


K 10 


150,000 


0 

max 


35 degrees 


K 11 


3,000 


0 

max 


10 deg/sec. 


K 12 


800,000 


8, time scaling factor 


1 
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Table II 



Potentiometer Expressions and Settings for Baseline Simulation 

(See Fig. 4) 



Potentiometer 



No. 


Expression 


Setting 


1 


e 0 Ki 

100K 2*«ax 


.9546 


2 


Vi 

100K 2 x max 


.9546 


3 


Simulation 

Technique 


.1 


4 


^4 x max 

l® z max3 


.1500 


5 


% z max 

100K-X 

2 max 


.1136 


6 


db 

4 


.05 


7 


db 


.05 


8 


multiple flow gain 


.5 


9 


multiple flow gain 


.5 


10 


Qmax (1st footnote) 

10 k 6 p l b 


.406 


11 


AL ®max^/ 180 ^ 2nd footnote ) 
20 % P Lmax 8 


.4654 


12 


P Lmax/p 

' r max 


.5000 


13 


^/Pmax 


.5 


14 


^Lmax 

K 12>Vax 


.0177 


IS 


L 

^max 


.8824 



26 



Potentiometer 

No. 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 



Expression 



Setting 



AP 



Lmax 



3.45X10 5 

^ll^max 

3.45X10 5 

^9 e max 

1.075X106 

W 



C3rd footnote) 



(3rd footnote) 



(3rd footnote) 



1.075X10 6 

( 3 - 45 X 105 / 1 . 12 9X106 ) 2 

(1.075X10 6 / 1>12 9 X10 6) 2 



(3rd footnote) 

(3rd footnote) 
(3d footnote) 



.6957 

.3044 

.9768 

.0233 

.6934 

.9066 



KlO .1286 

T 

Fmax 

4 

Cd R (1 . 129X10 6 ) (3rd footnote) .8714 
^Fmax 



.0010 



TFmax 

^max 1 



.1282 



26 

27 

28 
29 



Tpmav (Tc/Ts) 

iooe m is 

max 

n 

max 

9 max 

®max (4th footnote) 

Too 

Kg 

10 



.1269 



.2857 



.1750 

.25 
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Potentiometer 

No. 



Setting 



Expression 

30 ^7 .15 

10 

31 T Hmax .3370 

lOTFmax 



32 



Tpmax 

"^Frnax 



.6175 



1. This expression was later found to be incorrect. See the text. 

Initial Investigation. 

2. This expression was later found to be incorrect. See the text. 

Initial Investigation. 

3. Numbers appearing in the numerators or denominators of these expressions 
represent scaling values calculated by the NSRDC. 

4. Adjustment for SIN - COS generator design requirements. 
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these parameters modified the fundamental frequency of the oscillation, 
but did not eliminate it. It was found that the oscillation could be 
eliminated by the artificial introduction of viscous damping in the 
system, but only in much greater amounts than could be justified in view 
of apparent levels of viscous damping present in the actual system. No 
provision for viscous damping was included in the friction torque equation, 
since no source was immediately apparent in the structure, 
b. Suspected Origins of the Oscillations 

After some study of the baseline simulation, the NSRDC felt 
that the problem was located in either their recording equipment associated 
with their analog computer, in the mathematical models which had been 
developed, in settings of the simulator parameters, or in the basic assump- 
tions used in the derivations of the governing equations, such as a lack 
of understanding of the available sources of damping in a submarine 

4 

steering control system. 
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II. INITIAL INVESTIGATION 



The initial investigation was conducted in order to gain a more 
complete understanding of the simulation, and to determine problem areas 
requiring extensive investigation. The baseline simulation was duplicated 
on a large analog computer, techniques and models currently used in simulat- 
ing hydraulic control systems were studied, and information was gathered 
from personnel familiar with the baseline simulation development and sub- 
marine steering control systems. The baseline simulation was continually 
evaluated in terms of the information which had been gained, and several 
problem areas were found which required further analysis. 

A. DUPLICATION OF THE BASELINE SIMULATION 

The baseline simulation required duplication on local equipment in 

4 

order to provide a check of the correctness of the programming completed 
at the NSRDC. The Naval Postgraduate School's COMCOR CI-5000 analog 
computer was used in this effort. The CI-5000 is a general purpose 
computing system, including analog, digital control, and patchable logic 
sections, and incorporates enough components of various types to accommo- 
date the baseline simulation. 

The writer first attempted to directly adapt the baseline simulation 
analog schematic directly to the CI-5000, adjusting the simulation only 
to compensate for CI-5000 peculiarities, thereby assuming that the equations 
and potentiometer expressions of the baseline simulation were correct, and 
that the analog schematic diagram accurately modeled the governing equations. 
This first attempt was unsuccessful; analog components rapidly overloaded, 
and the system failed to reproduce the sample curves that the baseline 
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simulation at the NSRDC had produced when operated under identical 
conditions. It appeared that the baseline simulation's governing 
equations, potentiometer expressions, and analog programming required 
detailed checking for mathematical and functional correctness. 

The governing equations were found to be consistent with the sign 
conventions evidently adopted by the developer (See Figure 1) , with the 
exception of equation (7) , 



Q - -fv + k 6 p l • 

Examination of the geometry of the system (Figure 1) revealed that the 
sign of the y - term should have been positive, since a positive flow 
rate would cause a positive change in y, if the pressure drop were 

held constant. Also, for a zero flow condition (Q = 0) a positive change 
in y would result in a negative change in the pressure drop. Thus, 

4 

equation (7) should have read 

Q = Ay + K 6 P l . (7a) 

The NSRDC later confirmed that this was the case. 

Potentiometer expressions were found to be consistent with the baseline 
simulation schematic diagram and the governing equations, with the exception 
of potentiometers 10 and 11, which were associated with integrator 9 of 
the baseline simulation (Figure 4) . A term in the denominator of each 
potentiometer expression apparently had been dropped; the denominators 
should have read 10 (Pyj^) [K^+ (A^/K^) ] instead of lOfK^Py,^) and 

20(KgPj jIIiax ) for pots 10 and 11, respectively. The need for the additional 
term in each expressions denominator was later confirmed by the NSRDC. The 
additional terms modified the values of the two potentiometers significantly; 
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the value of pot 10 changed from 0.6197 to 0.4060; the value of pot 11 
changed from 0.4654 with a gain of 20 on integrator 9 to 0.6060 with a 
gain of 10. This resulted in an effective reduction of the system hydrau- 
lic compliance and in a change of the overall characteristics of the 
pressure torque. Tests later showed that the omission of these coefficient 
terms was the primary reason for the initial failure to duplicate the 
baseline simulation. 

Except for the two potentiometer settings, the analog program (Fig. 4) 
appeared to accurately represent the verified governing equations . Accord- 
ingly, another attempt was made to adapt the program to the CI-5000. Only 
slight deviations from the baseline schematic diagram were necessary in 
order to accommodate the CI-5000 component peculiarities. 

When the adapted baseline simulation was operated under the potentio- 
meter conditions listed in Table II, with the exception of pots 9 and 10 

4 

which were set as previously modified, the system behaved in the same 
manner as it had during operation at the NSRDC. Oscillations identical 
to baseline forms were observed in 6 and (Figure 5) . As a result of 
this duplication of performance, using several different types of recording 
equipment that exhibited varying response characteristics, it appeared that 
the recording equipment at the NSRDC was not the source of the oscillations. 
Further, since all errors in mathematics and in the settings of the simu- 
lator parameters required correction before the baseline performance could 
be successfully duplicated, it appeared that the errors which were found 
were the result of faulty communication of the baseline simulation to the 
writer, and were not present in the program used to generate the sample 
curves. Apparently, then, the mathematical models were consistent in 
sign conventions and had been manipulated correctly in the development 
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FIGURE 5 

BASELINE SIMULATION RESPONSE 
© 0 = 35 DEGREES 
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of the analog board, the potentiometers were correctly scaled, and the 
analog wiring was correct. Thus, the simulation appeared to be correct 
within the framework of the fundamental assumptions involved in the mathe- 
matical formulations; however, the oscillations were still present in the 
output . 

Two situations seemed to be possible. First, perhaps the oscillations 
evidenced in the NSRDC data were not representative of the behavior present 
in real submarine steering systems, and the simulation was therefore in- 
adequate. This situation could have occurred because fundamental errors 
or omissions were made during the formulation of the baseline simulation 
governing equations , say by the use of faulty assumptions regarding the 
significance of one or more physical processes present in a steering 
system. The situation could also have occurred because system materials, 
dimensions, and properties were improperly represented in the development 

4 

of the numerical values used in the "typical" system simulation, by the 
inclusion of inconsistent units conversions, unrealistic dimensions, or 
faulty physical materials constants. 

Second, perhaps the oscillations appearing in the NSRDC data were 
representative of real system behavior, and the difficulty had its origins 
in the harboring of an incorrect knowledge of the actual performance of a 
ship's steering control system. No actual system recordings were sub- 
mitted with the baseline simulation, and therefore no immediate comparison 
was possible between the simulation performance and real conditions in the 
system. 

B. BACKGROUND INVESTIGATION 

1 . Study of Currently Used Models 
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Assuming that the oscillations observed in the baseline simulation 
were not representative of real steering system behavior led to a back- 
ground investigation of techniques and models currently used in simulating 
hydraulic systems in order to develop a clear knowledge of the fundamental 
assumptions made during the development of the baseline simulation. This 
background examination consisted of studies in the hydraulic controls 
field, and communication with personnel familiar with the simulation and 
ship steering control systems. 

The search for currently used models and techniques produced several 
results. It was found that servovalves were usually designed so that they 
were the controlling element in hydraulic systems. However, in the base- 
line simulation, any error signal representing more than four degrees 
resulted in saturation of the control valve, resulting in the servovalve 
having no variable control over the velocity of the load, until the error 

4 

signal again fell below a value representing four degrees. Thus, the 
servovalve was virtually being used as an open-shut valve much of the 
time. The presence of deadband and multiple flow gain in the valve 
tended to reduce the sensitivity of the valve to small errors, and to 
increase it for errors in excess of four degrees, thus heightening the 
effect. 

It was also noted that if the expressions for the baseline simula- 
tion were transformed into Laplace notation, and deadband and multiple 
flow gain effects were temporarily neglected, the resultant expression 
would be 

Z (s) = ( K V K - 2 — \ E(s) (14) 

\ s+K 3 /K 2 / 
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FIGURE 6 

BLOCK DIAGRAM OF A TWO -STAGE SERVOVALVE 

WITH FORCE FEEDBACK 



where Z and E represent the Laplace-domain variables for z, the valve 
spool position, and e, the error signal, respectively, and s is the 
Laplace transform variable. This function was identical to one developed 
by Morse [1] for a valve incorporating centering springs instead- of force 
feedback, for use in systems where the load was assumed to be very small 
compared with the available hydraulic power. The actual valve being used, 
however, consisted of a two-stage, force-feedback design being used to 
control a load of the same order of magnitude as the available hydraulic 
power [2] . Merritt [3] developed a transfer function block diagram for 
a two-stage, force-feedback servovalve, which incorporated the general 
form shown in figure 6. A vast degree of simplification of Merritt's 
expression obviously would have been necessary in order to arrive at the 
expression used in the baseline simulation. Quite possibly, the co- 
efficients indicated by Merritt in many instances were of such magnitudes 

4 

that their effects were not important in the actual performance of sub- 
marine steering system servovalves, since high frequency effects would 
become unimportant when other system frequencies were quite low. Further- 
more, the pressure feedback to the flapper valve that Merritt included in 
his expression was possibly unimportant in this particular case. However, 
it seemed advisable to check the numerical value of each term in the 
detailed expression, in order to evaluate the validity of the baseline 
simulation’s expression. 

Through communication with Mr. Guy Johnson [4], the developer of the 
baseline simulation, it was learned that many of the servovalve design 
values, needed to calculate the coefficients in the Merritt block diagram, 
were difficult to obtain. It was also found that the expression used in 
the baseline simulation had been developed by a combination of experimental 
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and intuitive means, and that the servovalve simulation was considered 
to be sufficiently accurate for the particular design studies in which 
it would be used. Finally, it was determined through observations of 
the performance of the adapted simulation program that the oscillations 
did not appear to be affected by the servovalve action. Thus, further 
investigation of the validity of the servovalve simulation did not appear 
to be justified until the origin of the oscillations had been fully 
investigated. 




Figure 7. Functional Block Diagram of a Submarine Steering System 
Figure 7 is a functional block diagram of a submarine steering 
control system. The control system was not entirely closed loop in 
construction, since the error signal controlling the function of the 
servovalve was generated by the ram position rather than by the actual 
rudder position. This placement of the feedback element was consistent 
with established design practice [5] , since the inclusion of the compliant 
structure of the rudder and linkage within the servo loop could result in 
significant instability problems. Additionally, an electrical feedback 
element located outside the pressure hull of a submarine, in a location 
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FIGURE 8 

CURVES OF SEVERAL PARAMETERS 
OF Cl -5000 SIMULATION 
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FIGURE 9 

"frictionless" LOAD CONDITION,! =0. 
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flooded with sea water, would pose several design problems, including the 
need for increased reliability due to its inaccessability when at sea. 
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Figure 10. Circuit Modification- of Fig. 4, 
Introduction of Viscous Damping 
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(Figure 8) generated by the CI-5000 simulation revealed that the oscilla- 
tions were present only in the linkage and the load portions of the simulator. 
Manipulations of the baseline simulation friction coefficients failed to 
eliminate or modify the effect, although total exclusion of any friction 
in the system resulted in very small undamped oscillations in 0, and appre- 
ciable oscillations in 6 and after the control surface had essentially 
reached the desired position (see Figure 9). Further tests showed that 
small values of negative feedback across integrator 25, which were intro- 
duced through the use of the circuit modification indicated in Figure 10, 
resulted in elimination of the oscillations. This feedback across the 
integrator would be representative of damping in an actual system con- 
taining such elements as dashpots. With no other friction in the system, 
a feedback level of 0.005 ® was sufficient to eliminate the limit cycle. 
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FIGURE 12 

VISCOUS DAMPING LEVEL .12a ,T p =0. 

°max h 
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Furthermore, feedback levels of 0.12 2 — were sufficient to approach a 

9m 

critically damped condition in the load and linkage simulator. (See 
Figures 11 and 12.) If friction levels present in the baseline simula- 
tion were programmed into the CI-5000, approximately the same value of 6 
-damping as in the frictionless configuration was required to produce 
critically damped behavior, while the sustained limit cycle was eventually 
eliminated by the presence of the baseline friction itself. 

Damping terms proportional to 6 were totally absent from the base- 
line simulation equations. Intuitively, however, some damping would be 
present in any system in which viscous fluids such as sea water are dis- 
placed by moving objects, since work must be expended to keep objects 
moving at a constant rate, whether linearly or rotationally, through a 
real fluid. Thus, the neglect of damping in the linkage and rudder portion 
of the simulation did not appear to be realistic. ,A superficial examina- 
tion of the system did not reveal any large sources of damping, or even 
sources which would provide a fraction of that needed to approach a criti- 
cally damped condition in the linkage. Anticipating that a lengthy 
investigation of system damping would be necessary, the writer identified 
this area as requiring first priority for further investigation, and post- 
poned further analysis of it. 

Correspondence with Mr. Johnson resulted in the writer receiving 
some drawings of a 671-class submarine's steering control system [7]. 
Examination of these drawings revealed that functional concepts of several 
features of the design had undergone considerable alterations in the pro- 
cess of the development of the idealized model. A compensated actuator 
had been assumed in the simulation development; i.e., equal areas had 
been assumed on both sides of the actuator ram (Figure 1) . Examination 
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of the drawings, however, and discussion with Mr. Johnson revealed that 
the actuator on a typical submarine was in reality uncompensated; i.e., 
unequal areas were available on opposite sides of the actuating ram 
piston. The writer was informed, however, that the difference in per- 
formance parameters due to the neglect of the unequal areas was not of 
major concern in the overall reasons for simulating the steering system. 

Further study of the drawings revealed that the crosshead configu- 
ration had been simplified in such a manner that it was no longer represen- 
tative of the actual design. The crosshead, originally assumed to operate 
in a sleeve open at both ends, Figure 1, was in reality supported by a 
sleeve open at only one end in most submarine designs. Figure 13. The 
effect of this variation in the crosshead configuration could only be 
determined through a detailed knowledge of the dimensions, design, and 
clearances involved in an actual crosshead, and these parameters could 

4 

not be gleaned from the drawings that the writer had available. Further 
investigation of the crosshead influence on the system seemed to be justi- 
fied, since the assembly resembled a dashpot in appearance, and possibly 
was capable of providing the source of damping which seemed to be needed 
in the system. Consequently, a detailed analysis of the effect of a . 
typical crosshead structure on steering system performance seemed to be 
in order, and was identified as an area for further work. 

2. Verification of Typical System Values 

Several nondimensional coefficients, load parameters, and properties 
had been used to develop the simulation, such as the hydraulic resilience 
coefficient, the rudder mass, and inherent materials properties. Inaccurate 
development of system numerical values could have contributed the system 
oscillations. Attempts were made to determine the values of many of these 
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CROSSHEAD 




46 



FIGURE 13 

TYPICAL CROSSHEAD DESIGN 



constants from basic principles, but the effort was thwarted due to a lack 
of information regarding dimensions, physical constants, system construction 
details, and algebraic expressions for any of the nondimensional constants 
Kj through t> ecame apparent that much time could be spent trying to 

determine the validity of each derived constant without having the appro- 
priate algebraic expression available for examination. Therefore, this 
possible source of inaccurate system simulation was identified for further 
investigation, and was set aside. 

3. Actual System Performance 

Actual system recordings of 0, 0, and P were mentioned by Mr. 

J-i 

Johnson as not demonstrating the oscillations found in the simulation. 

The writer requested that such recordings be forwarded to him in order to 
gain a first-hand knowledge of the actual system performance, but was 
unsuccessful in obtaining any such recordings. Without any more information 
about the actual system performance than the available qualitative descrip- 
tion of the behavior, no reasonable conclusions could be reached about the 
actual system performance. Because an accurate representation of a system 
could not be claimed until a favorable comparison with accurate actual 
system performance could be accomplished, it was felt that such a determina- 
tion of steering system performance would be beneficial, and the subject 
was identified as requiring further investigation. 

C. SUMMARY OF THE INITIAL INVESTIGATION 
1 . Conclusions 

As a result of the duplication of the baseline simulation, and 
subsequent experiments performed with the program, the following conclu- 
sions were reached: 
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(1) The recording equipment at the NSRDC had accurately 
recorded the performance of the simulation, and was not the source of the 
oscillations . 



(2) The equations and potentiometer settings were consistent 
within the framework of the fundamental assumptions involved in the 
simulation's mathematical development. 

(3) The oscillations were confined to the load portion of 
the simulator, and did not appear to be sustained by overall system 
instability. 



n 

(4) A negative feedback level of 0.005 t— was sufficient to 



e 



m 



eliminate the limit cycle oscillation, and 0.12 -r— achieved a critically 

®m 

damped condition, in the load simulation under otherwise frictionless 

conditions. A negative feedback level of 0.12 was also sufficient to 

^m 

create a critically damped condition in the load simulator under baseline 

« 

friction conditions. 

Studies, and communications with personnel familiar with the simula- 
tion and with actual submarine steering control systems, led to the 
following conclusions: 

(1) The servovalve equation used in the baseline simulation 
was not developed to completely model either the type of servovalve used 
in the real system under study, or the type of load which was present in 
the actual system. The net effect of the substitution of a simple ex- 
pression for a complicated one was not immediately determined, except 
that it did not appear to affect the system load oscillations. 

(2) The actuator had been modeled in the system simulation 
as a compensated ram, rather than being modeled as the uncompensated 
actuator used in real systems. Although the effect of this simplification 
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on the system was unknown, communication with experienced personnel seemed 
to indicate that the effect was unimportant. 

(3) The actual crosshead design had been modeled in such a 
way that damping effects possibly present in the component were not taken 
into account. Further study seemed warranted. 

(4) Viscous damping effects in the load had been totally ignored 
in the development of the simulation, although such effects controlled the 
unwanted system oscillations, and normally would be present to some degree 

in any system incorporating the movement of structures through real fluids. 

A search for viscous damping sources in the load was assigned first priority 
for extensive investigation. 

(5) The origins of many of the nondimensional constants used 
in the simulation were obscure, and efforts to verify their correctness 
were unsuccessful, due to the lack of appropriate algebraic expressions, 

4 

dimensions, and materials constants. 

(6) Actual documentation of real system performance in the 
situations being modeled appeared to be vague and incomplete. Little 
information was available concerning the actual response characteristics 
of submarine steering control systems to step inputs. An extensive inves- 
tigation of real system performance appeared to be necessary. 



49 



III. LOAD AND LINKAGE DAMPING 



Of the areas identified in the initial investigation as requiring 
further research, only the search for sources of linkage and load damping 
was pursued in the present investigation. Analysis of the submarine 
system identified three locations where damping might possibly be present: 
the system's dynamic seals and bearings, the crosshead structure, or the 
rudder itself. No other components appeared to be designed in such a 
manner that they could possibly produce an appreciable quantity of damping. 

In order to determine levels of damping in the various components, 
detailed information was needed concerning dimensions, clearances, materials, 
constructions, and lubrication for the linkage parts and rudder. The writer 
was referred [8] to Mr. Leeland Smith [9] of the Mare Island Naval Shipyard, 

4 

who provided detailed drawings of the steering control system for the SSN 
671 class submarine. 

A. BEARINGS AND DYNAMIC SEALS 

The possibility was considered that the expressions employed for the 
friction present in the bearings and dynamic seals (equations 9 and 10) 
were not representative of the actual processes present in the submarine 
system, and that the inherent friction of the bearings and seals was 
better described as displaying viscous behavior. 

Keller [10] indicated that compression packings, such as those used 
at the pressure hull penetration and the actuator-piston rod interface, 
are normally lubricated by allowing very slight amounts of leakage past the 
seal. Examination of the SSN-671 drawings revealed that the pressure hull 
packing was additionally lubricated with grease. Despite these indications 
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that viscous effects could be present in these seals, discussions with 
Mr. Smith, William Kendall [12], and Leonard Chase [12] of the Mare Island 
Naval Shipyard resulted in the writer's learning that friction in such 
seals was best approximated as sliding friction, since friction loads were 
normally observed to be fairly independent of actuator rod velocity, due 
to the velocities and allowable leakage rates present in the system. Thus, 
the baseline model of seal friction appeared to be justified. 

Study of the rudder bearings seemed to support the use of the baseline 
simulation's bearing friction expressions. The large rudder weight, com- 
bined with hydrodynamic loads, small rates of rotation present in the 
rudder stock, and the ill-defined lubricating qualities of sea water, 
resulted in enough analytical uncertainties being present to discourage 
any attempt to model such bearings as either slipper-type or journal-type 
surfaces. Again, the baseline simulation expression seemed to provide a 
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bearing surface description that probably was as accurate in magnitude and 
parameter dependence as any other expression that could be analytically 
developed using known conditions and existing lubrication theory. Experi- 
mental and analytical work of a scale beyond the scope of this project 
could prove useful in future design work involving bearings and seals of 
the types present in the submarine system. 

B. CROSSHEAD 

1. General Description 

Examination of the SSN-671 drawings, and discussions with Mssrs. 
Smith, Kendall, and Chase concerning general construction details of 
crossheads, revealed that, as a rule, crossheads somewhat resembled dash- 
pots in construction, due to their piston-in sleeve arrangements, and their 
free-flooded nature. Provisions were normally made in the construction of 
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crossheads to eliminate any possible damping effect by the drilling of 
large holes through the piston portion of the structure. Because damping 
sources were being sought, the writer decided to analyze the crosshead in 
order to determine what damping effects, if any, were present. 

Figure 13 shows a crosshead design based on information taken from 
the SSN-671 drawings. The crosshead consists of a partially bored out 
cylinder 19 inches in length, with a diameter of 11.986 inches, and four 
one-inch I.D., eight-inch in length holes axially bored through the solid 
head of the piston. This piston rides in a twelve-inch I.D. tube, typically 
sealed along its length except at the rear end, which is open to allow the 
passage and free rotation of the link connecting the actuator shaft ex- 
tension and the rudder tiller. The piston incorporates several lock- 
preventing 0.75 in. X 0.125 in. circumferential grooves, and two 0.125 in. 

X 0.5 in. axial slots along the outside radius. As a result of the radial 
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clearance between the piston and the sleeve, the slots, and the holes in 
the crosshead, the entire sleeve is free flooded, and flow of sea water 
past the crosshead must take place as it is moved in and out of the sleeve. 

2. Assumptions 

In approximating the behavior of the flow past the crosshead and 
the viscous drag associated with that flow, it was assumed that the fluid 
velocities through the holes, slots, and radial clearance of the crosshead 
were considerably greater than the motion of the piston relative to the 
sleeve, thereby implying the existence of fixed orifices of irregular 
shape. Laminar flow was also assumed. 

3. Equations 

a. Eccentric Piston in a Sleeve 
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For an annulus formed by a solid eccentric cylinder in a sleeve 



with no relative movement involved, the flow rate can be expressed as [13] 



irrc 
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' -mv 



( p r p 2) 



(15) 



where r is the radius of the sleeve, y is the absolute viscosity of the 
fluid, 1 is the length of the tube, e is the eccentricity of the piston with 
respect to the center of the sleeve, and (Pj-P2) is the pressure drop across 
the length of the tube. 

b. Laminar Flow Through Rectangular Passages 
For the slots [13] , 
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(16) 



where w is the width of the slot, h is the depth, and all other symbols 
retain their previous meanings. * 

c. Flow Through Short-Tube Orifices 

Merritt [14], combining the work of Langhaar [15], Kreith [16], 
and Shapiro [17] , developed an expression for the flow rate Qq through short 
tube orifices: 




(17) 



(18) 



where C is the loss coefficient, p is the fluid density, D is the hole 
diameter. Re is the Reynolds Number based on the mean velocity of the 
flow in the tube and the hole diameter, and all other variables retain 
their previous meanings. If equations 15, 16, and 17 are manipulated 
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briefly, the pressure drop can be expressed as a function of the flow 

rate Qq, the kinematic viscosity v, D, and p: 



P l- P 2 = 



8p 
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If a total mass balance is made across the crosshead, it is found that 



Ach y = MQq + NQ S + Q c (20) 

where A^ is the cross-sectional area of the crosshead, M is the number of 
holes in the crosshead, and N is the number of slots. 

The force caused by the pressure difference across the crosshead 
can be expressed as an equivalent torque at the rudder pivot, by multi- 
plying the crosshead force F = \;H^ P 1 -P 2 ^ by the length L of the tiller 
arm: 

T CH - \h l tP r P 2) ' < 21 > 

Figure 14 shows curves of versus y generated through the solution of 
equations (18) -(20) using the previously listed crosshead dimensions, p 

3 r o 

equal to 64 lbm/ft , and v equal to 1.41 X 10 -D ft /sec (water at 50 degrees 
F) , for various sizes of piston holes, assuming four holes and two slots, 
and an eccentric piston in the sleeve. To generate the curves, equations 
(18) and (19) were solved for as a function of Qq. Then Qg and Q^. 

were calculated, and appropriate multiples of Qq, Qg, and Q^, were used in 
accordance with equation (20) to produce y. Equation (21) was then used to 
determine T CH 

For any crosshead incorporating four holes any larger than 0.1 
inches in diameter, it was found that zero pressure differential existed. 

As the diameter of the holes was decreased, with the dimensions and number 
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FIGURE 14 

CROSSHEAD TORQUE 
FOR VARIOUS HOLE SIZES 
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of slots and the parameters for the annulus remaining constant, a pressure 
differential began to exist across the structure. Finally, as the holes 
were eliminated, (P 2 /P 2 ) approached levels sufficient to create significant 
values of Tqj. For the values of the zero-holes curve used in Figure 14, 
T CH could be expressed as 



T CH = 7 - 48X1 ° 3 y 



y e 

or, since - — >\> -> — , 
^m ®max 



T C H = (7.48X10*ymax) C^) 



( 22 ) 



(23) 



If the above equations were scaled for introduction into the 
simulation, potentiometer 33 would be set according to the expression 



7. 48X10^)011 ax 



(24) 



33 100 I 0 max 

For I = 9100, and y max = 3 inches per second, 4 

a 33 = .00412. 

Since this value of the potentiometer represented the maximum available 
Tqj, it was seen that the actual effect of the crosshead on the system 
damping was very small. In fact, using one-inch diameter holes would make 
the crosshead damping essentially zero. However, since a feedback level of 

Q 

.005 -» — was previously found to be sufficient to eliminate limit cycles in 
e m 

0 

the "frictionless" system simulation, .00412 could assist other small 
sources in the damping of the oscillations in the baseline simulation. Al- 
though the results of the investigation of the crosshead damping were less 
spectacular than had been desired originally, since the 0- feedback sufficient 
to create a critically damped condition in the load had not been justified, 
the author felt that the identification of a small source of damping had 
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been achieved, that identification possibly proving to be of some worth in 
future design work if a small amount of damping somehow proved necessary. 

Some argument could have been made that further modification 
could create more damping. However, for this to be done, the system would 
need to be sealed at both ends, due to the approach of ?2 to the vapor 
pressure of water in the contained volume when the crosshead was moving 
out of the tube. Of course, the open end of the crosshead sleeve was not 
capable of being closed, due to the primary function of the crosshead, 
which was to provide a portion of a slider-crank mechanism. The previously 
developed levels of Tq^ thus appeared to be the maximum practicable. 
Additional damping could be generated, if needed, through the use of a 
large dashpot in parallel with the actuating linkage. However, such a 
dashpot was not apparently needed, as was shown through investigation of 
the inherent damping of the rudder. 



C. RUDDER 

The description of the moments involved in the rotation of a rudder 
about its pivot, while moving through a fluid, would in general include 
functions involving the zero, first, and second-order derivatives of rudder 
angle with respect to time. That is, if a moment balance were written for 
a rotating rudder, the equation could in general take the following form: 



M = fjCtf) + f 2 (0) + f 3 (0) 



(25) 



where f , f ? , and f would be functions of some yet-to-be-determined form. 

Further consideration of the natures of f ^ and f^ revealed that expressions 

for these two functions were present in the baseline simulation. If f 2 = 

f = 0, the general expression is reduced to 
o 

M = f x (V ) (26) 
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( 27 ) 



and is recognized as being of the form 




where I is the polar moment of inertia of the rudder and its fluid added 
mass. If f^=f2=0, f^ is recognized as the sum of all 0 - dependent 
parameters in the rudder, including the hydrodynamic torque caused by the 
rudder lift and the difference between the location of the rudder's pivot 
and its center of pressure, and the 0 - dependent portion of the system 
friction. 

The function f^Ce) was assumed equal to zero in the baseline simulation; 
that is, no 0 - dependent effects were included in the rudder's equation of 
motion (Eqn. 8). The validity of this assumption was examined in the final 
portion of the investigation of load and linkage damping. 

1. Existence of a Non-Zero Function f 2 (® ) 

Initial considerations of the problem involved the search for any 
indication that a non-zero function f 2 (© ) even existed. Intuitively, such 
a relation should exist in some magnitude, at least under a condition where 
the rudder was subjected to a constant turning rate while immersed in a 
stagnant, viscous fluid, since power would necessarily be expended to keep 
such a structure moving under these conditions, thus resulting in some level 
of torque being applied to the rudder shaft. This applied torque, then, 
would by Eqn. (25) be equal to f 2 (9 ) » and therefore f2 could not be identically 
zero. 

Another indication of the existence of a non-zero f2 was provided 
by Landweber [18] , who indicated that forces on lifting surfaces in unsteady 
motion were due to a combination of inertia effects and circulation, and 
that a transient variation did exist in the circulation around a thin air- 
foil impulsively given an angle of attack a, causing a transient behavior 
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of the lift coefficient of such a surface. Additionally, Theodorsen [19] 
studied the forces on a foil undergoing harmonic oscillations in connection 
with studies of control surface flutter and aerodynamic instability, and 
found that the actual lift coefficient was related to a quasi-steady lift 
coefficient corresponding to the instantaneous angle of attack of the foil, 
and time. Thus, the existence of a non-zero function f 2 (e), while not 
confirmed by the above references, was at least shown to be plausible. The 
writer therefore attempted to locate any work that had been done in the 
analysis of the forces and turning moments on airfoils or rudders subjected 
to constant turning rate. 

. 2. Damping in a Rudder Immersed in a Stationary Fluid 

No references were found describing the resisting torque developed 
by a rudder subjected to a constant rate of rotation while immersed in a 
viscous fluid. Accordingly, the writer attempted to approximate the drag 

4 

which could be expected, using elementary concepts and simplifying 




Figure 15. Rudder Construction Assumed For 

Stationary Rudder Torque Calculation 
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The rudder was approximated by a series of very small flat plates, 
assumed to be separated from each other so that the flow could move around 
the plates (Figure 15) . For a flat plate placed in a stream flowing perpen- 
dicular to the plane of the plate, the drag force F^ can be expressed as 




(28) 



where in this case is the drag coefficient, p is the density, is the 
area of the small plate, and the velocity V is related to the radius r and 
the angular rate of rotation 0 by V = r9. Thus, for a differential area 



dAp, 




pr 2 0 2 



dA p 



If the area is assumed to be equal to the span of the rudder S times the 
incremental distance dr,dA^ = Sdr, then 



dF n = C r 



( 



psV \ 



Sdr 



The incremental torque dT^ caused by developed at a distance r from the 
pivot point of the rudder would be 

dl R - C D ( " 4 - ) Sr3dr (29) 

Thus, for the portion of the plate from r = 0 to r = r , 



x C ° PSr ° 4 
T R = —8— 6 



(30) 



For a flat plate in a normal stream with a large S/dr ratio, is 
approximately 1.95 [20]. For the two rudders of a typical submarine, 
the following physical parameters are representative: 
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Parameter 



Value 



S 

r , rear portion 
r Q , forward portion 
P, sea water 



10 Feet 
7 Feet 
3 Feet 
64 lbm/ft 3 



The entire 6 - dependent torque would be equal to the Siam of the effects 
of both upper and lower rudders, including forward and aft portions: 



T e - 2 ( t Rforward + Tr aft } • ( 31) 

Substituting in the previously mentioned values, it is found that 

* O 

Tq = 88 e , in.-lb f (32) 

• • 

when 0 is expressed in degrees per second. If 0 = 10 deg/sec, T- = 8800 

uITlcLX 

in.-lbf. If an asymmetric square analog component and a potentiometer were 
placed in the simulator as in Figure 16, this source of damping could be 




Figure 16. Analog Circuit Diagram for Introduction of Tq 
T he appropriate potentiometer setting expression for simulating T^ would be 

T^m 
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(33) 



1000 18 
m 
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Substituting numerical values. 



a 34 = .968 X 10‘ 3 . (34) 

Thus, the initial estimate of available rudder drag at zero velocity indicated 
negligible amounts of damping were present in the rudder turned at a constant 
rate in a stationary fluid. The writer felt that the results were incon- 
clusive, however, since the assumptions which were made in the calculation 
of t r were very rough indeed. The rudder was not formed of many small strips 
with spaces between them; nor was the flow as simple as that described by the 
model. Such a complicated flow as would be present along the surface of a 
large flat plate rotating in a quiescent fluid required much more analysis 
and study than- the writer afforded it. He therefore concluded that future 
studies in this area were necessary, and set aside consideration of the 
stationary rudder case. 

4 

3 . Damping in a Rudder Immersed in a Moving Stream 

A literature search revealed that an analysis had been made of the 
torque behavior of a turning rudder immersed in a moving stream. Okada [21] 
published a paper dealing with the effect, which was subsequently expanded 
by Kennard and Leibowitz [22] in their studies of rudder performance as 
affected by ship turning, steady change of rudder angle, and propellor race. 

Kennard and Leibowitz [23] noted that the twisting moment observed 
on a rudder shaft while the rudder is being turned may differ considerably 
from its value when the rudder is stationary. They went on to explain that 
this effect was due to circulatory motions around the foil caused by 
viscosity effects of the water, which subsequently generate vortex sheets 
in the water downstream from the rudder's trailing edge, and cause moment- 
modifying variations in the pressure distribution over the rudder's surface. 
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Variations from steady state twisting moments by as much as thirty percent 
were mentioned as being observed in rudders turning at constant rates. 

The mathematical difficulties associated with the problem were 
simplified by several approximations. The rudder was replaced by the flat 
plate, shown in Figure 17, of uniform chord C and infinite vertical span. 
The water was assumed to approach the plate at a uniform rate and at an 
angle of attack a; its velocity was denoted by U. The rudder twisting 
moment per unit span was identified as M^, and was assumed positive when 
operating in a direction to decrease a. The pivot point of the plate was 
identified as H, and was located at a distance hC from the leading edge, 
h being assigned some numerical value between 0.25 and 0.50. The rudder 
turning rate was assumed positive for increasing angles of a. 




Figure 17. Flat Plate Approximation to a Rudder 
Turning at Constant Rate 

Mjj, the twisting moment per unit span, was placed in nondimens ional form 
by dividing it by (1/2) C^U^, to form the dimensionless moment coefficient 
Cj^. An expression for was then derived, with C being represented as 
the sum of a quasi-stationary part an <i a 0 - dependent part, AC^. 

Thus , 

Sm = C MH0 + AC MH t 35 -* 



63 



20 



30 



40 



50 




C< U 

r d<x 

dt 



FIGURE 18 
A VERSUS fm 



64 



where 



1 

Cj^q = (4 -h) 2 ira (36) 

ACmh = ( £(§ -h)(l-2h)-(I -h)Aj (37) 

and 

_ C da 
’[) dt 

and A is a function of Si and a . The effect of A is to increase the value 
of ACj^j as a increases, as is shown in Figure 18. 

Kennard and Leibowitz reported that the expression for was not 

entirely correct, due to the neglect of a nonlinear effect in its develop- 
ment. That nonlinear effect concerns the aft movement of the center of 
pressure in a rudder as a is increased, which causes a change in the sign 
of with increased a. Also, they indicated that values for AC^ obtained 

from model tests suggested that the calculated values of AC^ were somewhat 

* 

too high, although the scatter of the experimental data in the tests was 
considerable. They further concluded that the general trends resulting 
from turning real rudders in water were modeled in their theory, although 
values for AC^ were somewhat too high, such as a calculated increase in 
one case where h = 0.3 of roughly 40 percent, where in practice the increase 
due to AC^ ranged from 15 to 40 percent. It was suggested that tests with 
larger rudders would be worthwhile, and that refinement of the analysis was 
needed, expecially in including the effects of finite rudder spans. 

4. Application to the Simulation Problem 

The analysis just described seemed to be applicable to the simulation 
problem. Although the expression for was apparently unrealistic, the 

development ofAC^ was apparently useful in predicting the effect of rudder 
turning on the observed moments on a rudder, although calculated values 
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seemed to be somewhat higher than values which were experimentally obtained. 
In spite of its inherent shortcomings, incorporation of this modeling of 
the effect of rudder turning rate on the observed twisting moment seemed 
justified in light of the reported magnitudes of change in twisting moment 
observed in practice. 

In order to apply the above analysis to the simulation, several 
assumptions were made. It was assumed that a was identical to 9 ; that is, 
that the fluid flowed parallel to centerline of the submarine. Analysis of 
drawings of the SSN-671 submarine rudder revealed that h was approximately 
equal to 0.28; in order to compensate for the reported high nature of the 
calculated AC^, and to simplify calculations, h'was assumed to be 1/4, 
thus eliminating dependence on A. Also, the sign convention of was 
reversed to conform to the simulation's convention. The above assumptions 
when used in the expression for AC^ resulted in 

4 

Re-dimensionalizing to form AM^, the change in twisting moment per unit 
span due to 0, and then multiplying by the span length S, resulted in the 
expression 

AM h s = - J pC 3 SV0 . (39) 

For a SSN-671 rudder, the average value of C is about 10 feet, and S is 
also about 10 feet. In a typical submarine V can vary from zero to in 
excess of 35 ft/sec, and 0 was assumed in the baseline simulation to vary 
from zero to 10 degrees/sec. Assuming p = 64 lbm/ft 3 , the maximum value of 
the twisting moment, (AM^S) max based on Vo = 35 ft/sec, would be 

(AM H S) max = 5.72 X 10 5 in-lbf . (40) 
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Typically, a submarine has two rudders, of approximately the same design. 

The total twisting moment available would therefore be twice the above 

6 

value, or 1.144 X 10 in-lbf. An expression for aI%S/(AM^S) would be 



(“hSW 




(41) 



AMpjS/ (AM^S) max is dependent only on 0 for a given value of V, and therefore 
could be simulated as damping across integrator 25 as in Figure 10. The 
appropriate expression for the potentiometer setting 833 would be 



33 



C 4M H s WW v o) 



10 max (100)H 



(42) 



For V = V Q , and the conditions otherwise assumed above, since V appeared to 
be set equal to V Q for the other V - dependent parameters in the baseline 



simulation, 833 = 0.1259. A level of 0.12^- feedback was previously found 
to be sufficient to approach a critically damped behavior of the simulation, 
and therefore eliminate the oscillations for the system. Apparently, then, 
a theoretical justification for the inclusion of that much damping in the 
system had been found. The fact that AC^ was dependent on V did not allow 
the extension of the above damping levels to lower velocity, since as V 
decreases, so must the level of damping. Through communication with Mr. 
Johnson of the NSRDC, the author was led to believe that the damping effect 
somehow assumed a significantly large value even at low or zero velocities, 
thus indicating the existence of some unexamined or undiscovered contributing 
effects such as those previously discussed, which would be independent of V. 
More work in the area of rudder twisting moment and lift variance with rudder 
turning angle and other parameters obviously could be done. 
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5. Summary 

The investigation of linkage and load damping resulted in the 
location of a significant source of 9 - damping in the motion of a rudder 
immersed in a moving fluid. A significant amount of damping was evidenced 
by the large variation of the observed twisting moment produced by the 
turning of a rudder from that which would be observed if the rudder were 
held stationary at a given angle of attack. On the basis of rudimentary 
analysis, no large amounts of viscous damping were located in the packings, 
or the crosshead structure, although a slight amount of damping could be 
caused in the crosshead through redesign of the structure. A very rough 
estimate of the twisting moment observed on a rudder immersed in a stationary 
fluid indicated that little available damping existed unless the rudder was 
being moved through the fluid. Obviously, more analytical and experimental 
research needs to be done concerning the transient lift and twisting moment 

4 

characteristics of moving and stationary rudders in fluid flow fields. 
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IV. RECOMMENDATIONS FOR FURTHER WORK 



Several topics were identified as requiring investigation during the 
progress of this project. Due to time limitations, the writer was unable 
to pursue any of them in detail. Subjects uncovered in the initial investi- 
gation which remained significant at the end of the period that this report 
covers included the need for examination of the actual performance of sub- 
marine steering control systems under varying conditions, including a wide 
range of ship velocities. Studies of the effect of using a higher order 
servovalve and uncompensated actuator models in the simulation are necessary. 
Possibly, the re-checking of the consistency of the numerical values used 
in the baseline simulation is warranted, although the discovery of a signi- 
ficant source of system damping seemed to make this unnecessary. Additionally, 

4 

the investigation of linkage and load damping emphasized the need for con- 
tinuing research into unsteady rudder system performance characteristics, 
including rudder lift and twisting moment variations with such variables as 
rudder turning rate, ship velocity, ship turning rate, hull configuration, 
and propeller race. 
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V. CONCLUSIONS 



1. Analysis of unsteady rudder behavior identified a significant 
source of rate-dependent torque which had been neglected in the development 
of a submarine steering control system's analog simulation. Programming 

of representative levels of feedback in the simulator virtually eliminated 
oscillations in rudder rate and load pressure drop which were thought to 
be unrealistic. 

2. The actual performance characteristics of a submarine steering 
control system under varying operational conditions appeared to be loosely 
defined, and quantitative studies of system performance seemed to be 
necessary in order to continue the development of a realistic simulator. 

3. Studies showed that the inclusion of more complex models of the 

4 

servovalve and actuator could have important effects on the flexibility 
and realism of the simulator, and that further work should be done to 
accomplish this end. 

4. Continuing theoretical and experimental research in the field of 
rudder behavior under unsteady conditions was suggested in order to improve 
the current concepts employed in steering system design and simulation. 



70 



Appendix A 

Explanation of Notation 



1. Sat ( x ) uL 
LL 



.ul 



i 


ul 


x>ul 




x 


ll^x^ul 




) 11 


X<11 



This can be represented graphically as follows: 




Sgn(x) 



+1 


x>0 


0 


x=0 


-1 


x<0 
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